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We present an analytical formulation of the recent one-shot decoupling model of Bradler and 
Adami arXiv:1505.0284 and compute the resulting “Page Information” curves, for the reduced 
density matrices for the evaporating black hole internal degrees of freedom, and emitted Hawking 
radiation pairs entangled across the horizon. We argue that black hole evaporation/particle produc¬ 
tion has a very close analogy to the laboratory process of spontaneous parametric down conversion, 
when the pump is allowed to deplete. 


I. INTRODUCTION 


Recently there has been renewed interest in the use of a trilinear Hamiltonian to phenomenologically and explicitly 
unitarily model the evaporation of quantized internal degrees of freedom (treated as scalar bosons) of a black hole (BH) 
while commensurately generating Hawking radiation pairs one of which falls inward behind the BH horizon and 

the other, observed mode that has a thermal distribution, at least for early evolution times. The goal of representing 
the internal degrees of freedom of the BH as quantized is to be able to study both the evaporation of the BH as well 
as the effect on the thermality of the emerging Hawking particles. In the work of Nation and Blencowe [T] and Alsing 
[2] a single Hawking radiation pair mode was considered, that coupled to an internal quantized mode of the BH. Both 
models reprodnced information curves conjectured by Page HU] which predicted that the difference between the von 
Neumann entropy of an effective thermal density matrix for the outgoing Hawking radiation and the actual entropy, 
the so-called Page information I = S{pthermai) — S{p), would exhibit the following behavior, namely (i) for early 
evolution times I would be flat and nearly zero, since the BH is essentially un-depleted, and the emitted Hawking 
radiation has a thermal nature, and (ii) for long times, when the number of emitted Hawking particles is on the order 
of the remaining number of particle in the BH, I would begin to rise rapidly as the BH evaporates, and the emitted 
Hawking radiation deviates from a thermal state. 

In a recent paper, Bradler and Adami [3| have expanded upon the work Nation and Blencowe [1] and Alsing [2] 
and developed a one-shot decoupling model for BH evaporation, in which the quantized modes of the BH are emitted 
sequentially in time into the out-going and infalling radiation modes m- The authors use the method of lattice 
paths to compute the von Neumann entropy of the reduced density matrix for the BH and show that its resulting 
time evolution is again qualitatively similar to that predicted by Page, in which information emerges from the BH at 
a time when the BH has evaporated to roughly half its initial size (population). The new feature here is the use of a 
sequence of temporally emitted Hawking radiation modes (vs one) to more realistically model the creation of a train 
of Hawking pairs, each of which interacts with the BH quantized internal degree of freedom for some finite amount 
of time before the next subsequent Hawking pair emission event. Such a model conforms to the suspected physical 
generation process of Hawking radiation as suggested by many past HU] and recent works [5] . As discussed by Mathur 
[5] in his informative, pedagogical introduction to the BH Information Problem, Hawking pairs are created in a region 
of curvature distortion near the BH horizon, and subsequently propagate away from the region of generation, while 
the BH horizon shrinks in radius. As such, the next generated pair does not interact with the previously generated 
pairs, and is coherent only with the BH internal degrees of freedom for a time on the order of average time between 
emission events. This is what the one-shot decoupling model of Bradler and Adami H seeks to capture. 

The goal of this present paper is to analytically expand upon the one-shot decoupling model based on the trilinear 
Hamiltonian of Alsing H, and discussed in detail in the next section. The new aspect of this present work is the 
development of analytic approximations for the probabilities to the quantized BH / Hawking radiation state which are 
subsequently used to generate Page information curves, without the need for excessive computational resources. We 
further argue that the presented model of combined BH evaporation/Hawking pair production has very close analogies 
to the laboratory process of spontaneous parametric down conversion (SPDC) H [10] with a finite (vs infinite) BH 
‘pump’ source capable of depletion (the focus of H). 

The outline of this paper is a follows. In Section ]ll] we describe the trilinear Hamiltonian used to the manifestly, 
unitarily model the BH evaporation/Hawking pair production, and in Section HI review the one-shot decoupling 


model of Bradler and Adami |3]- Section IV we develop analytic approximations to the probability distributions 
for the reduced density matrix of the state of the BH, and for the state of the emitted Hawking radiation. Based 
on these approximations to the probabilities, we develop Page Information curves which illustrate the predictions 
conjectured by Page [5]. In the final discussion section Section|^we make the analogy of the this one-shot decoupling 
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to the process of SPDC generated by a laser powered by a battery of finite capacity. 


II. REVIEW OF TRILINEAR HAMILTONIAN MODEL 

In the work by Alsing [2] the evaporating BH hole was modeled by the manifestly unitary process given by the 
following trilinear Hamiltonian 


H - = T 


apalal 


d^Cls^i 


.)■ 


( 1 ) 


where ap are the creation and annihilation operators for the BH quantized, internal degrees of freedom, and 

aj. Os and at, are the creation and annihilation operators for the outgoing (observed) and in-falling (unobserved, 
behind the BH horizon) modes of the emitted Hawking radiation. Hp^-^ also models the process of parametric down 
conversion with a depleted pump source [111 I12j . Following the notation of [2] we use the subscript notation p to 
denote the BH ‘pump’ soruce, and the labels s and i to denote emitted the ‘signal’ and ‘idler’ modes respectively. 
Without loss of generality we take the signal mode s to be the particle and the idler mode i as the anti-particle 
emitted externally and internally, respectively to the BH horizon. Defining the initial state of the combined system 
as 


l'*/’(0)) ~ l''^po)p|0)s|0)j — |n.po)p|0)g_j. 


( 2 ) 


where ripo S> 1 is the initial boson occupation number for the BH. By solving the Schrodinger equation for the state 
IV’(t)) = J2n=o Cn('^) \n'po-n)p\n)s\n)i = X)n=o Cn('r) \npo-n)p\n),,- in both the short-time (upo » n) and long-time 
{npQ,n » 1) one obtains the solutions [2] 


= [(l-z)z”]^/\ 

Cniz,T) = [c<{f{z)/{l+f{z)),T)Y^^ 
' l + y/^ 

\fz 


f(z) = 4e“ 




0<z<z*, z* ^ 0.506407 for ripo > 1, 
z* < z < 1, 

z = tanh^(T), T = y/n^rt. 


(3) 

(4) 

(5) 


In the shorttime limit (ripo ^ n) the BH ‘pump’ source is essentially undepleted, and one can factor out the state 
using the approximation ripo — n « ripo, to yield the well studied two-mode squeezed state \(j))s^i form of the 
emitted Hawking radiation 


,r))< 

~ \f^po)p |</’)s,i) 

Z < Z*, 


(6) 

(sqzd) 


TIpO 

» Vl - z ^ z'^/‘^\n),-i, {upo oo). 

(7) 

s,i 

V I — z"pO + l 


where |n),y = (aj,a|)”/n! lO)^^. 

In [5] the initial number of bosons ripo in the BH was taken to be finite, though sufficiently large ripo 3> 1 so 
that without loss of generality all upper limits of summations could safely be taken as infinity (though in numerical 
computations, a finite value of ripo was utilized). In this work we will be especially careful to keep track of the upper 
limits of all summations, with UpQ large but finite so that Eq.(|^ is formally given by 


Kz,t) 


(1-^) 

(I - z"pO-l-l) 



VPn(z,T), 


( 8 ) 


such that 


n—O 

A physical motivation for utilizing the trilinear Hamiltonian Eq.Q to model BH evaporation/Hawking pair creation 
is as follows. The standard approach to modeling Hawking pair creation taken by Hawking and numerous subsequent 
authors is to treat the gravitation field of an eternal BH (typically taken as Schwarzschild, without loss of generality) 
as a classical background field, to which a quantized, scalar boson field is coupled [T3]. The subsequent state of emitted 
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Hawking radiation pairs was shown by Unruh HI to be a two mode squeezed state which upon tracing 

out of the (idler) mode i which falls behind the horizon, becomes a thermal state |15] . Very early on, researchers 
such as Boulware [5] and Gerlach [7], and more recently by authors such as Stojkovic et. al. HU and Alberghi et. al. 
HU, considered BH evaporation from a collapsing thin shell of matter. Here, the BH matter is still treated classically, 
and is characterized by a time dependent, shrinking horizon radius. The Hamiltonian derived for a scalar boson field 
coupled to this classical gravitation field (see HSl HZ] ) leads to a quantized harmonic oscillator with a (exponentially) 
time varying frequency. Such a Hamiltonian can be written in the form |18[ 119] H = l/2r^(r) (a^^ + a^) involving 
the generators of su(l, 1) and well known to generate single mode squeezed states [HI HH]- Here ^(r) represents the 
classical (i.e. c-number) ‘driving field’ of the collapsing shell of matter. The single mode squeezed state arises from 
the coupling of a single quantized scalar boson held to the classical gravitational held. If two bosons were coupled 
to the held, or a single, complex boson held, the Hamiltonian would be of the form H = l/2r^(r) (a^^ + as a^) 

where we have labeled the correlated emitted Hawking pair as a signal and idler mode. Again, the gravitation held 
shows up as a classical c-number driving held. In quantum optics, such semi-classical models are familiar, where 
the occupation number of a strong driving pump laser is so large, that for all intent and purposes concerning the 
subsequent particle statistics, it is eminently reasonable to consider it as a classical c-number held [9lHQ]. The next 
logical step, in order to incorporate the quantum statistics of the pump, is to replace the classical driving held ^(t) 
by quantized mode Up and a^p. This is the trilinear Hamiltonian Eq.(lJ) used in this work (and in HH3j)- Since at 
present, there does not exist a well accepted description of a quantizeclgravitational held (Schwarzschild, collapsing 
shell, etc...), this trilinear Hamiltonian, though physically motivated and eminently reasonable, must be taken to be 
at best phenomenological. It’s main advantage is that it manifestly unitary while capturing the essential qualitative 
features of generating across-the-horizon entangled Hawking radiation pairs, as well as serving as a simple model for 
BH evaporation. 


III. ONE-SHOT DECOUPLING MODEL 


The work of Bradler and Adami [3| generalizes the process used in Section |h] to the more physically relevant 
one-shot decoupling state in which over some period At a Hawking signal/idler pair is generated by the curvature 
in the vicinity of the instantaneous BH horizon (with radius formally proportional to n.p(T), the instantaneous BH 
occupation number), which then travels away from the region of generation, never to interact with the BH again, 
as discussed by [ZlIH] (and as suggested in the summary/conclusion of [5]). Such a model is physically analogous to 
spontaneous parametric down conversion (SPDC) in a nonlinear crystal of hnite length (as is typical in laboratory 
experiments) in which the entangled (here) photon pairs are generated inside the crystal (as some random spatial 
position). Upon exiting the crystal, the signal/idler pair no longer participate in the SPDC process described by 
the Hamiltonian Eq.Q. Thus the generated signal/idler modes emerge as a temporal sequence of emitted entangled 
pairs. Bradler and Adami model this using the initial state 


N 


1^(0)) = |npo)p n = Ko)p ® ® |0),^y^ 0 ... 0 |0),^y 


( 10 ) 




where r = NAr^ and N is the number of time slices. The evolution of the state |V'(0)) in Eq.(lO) is given by [3] 


N 


N 


|^(r)) = U(t, 0)|T(0)) |^'(0)) « e 




Wpo)p bl 


( 11 ) 


k=l 


k' = l 


where T is the time-ordered product and in the second equality we have used a simplified version of the Trotter 
expansion valid for N small time slices of size At, with Up^k = 6“*^'’'=^’’'=^'’’ acting on modes p and {sk,ik)- 
After the hrst time slice, the wave function is 


l'I^(l)) 


--pu , - iV 

ni—0 k'—2 

llpO , - ^ 

^pO -- 

|npo)p0 y |ni)i, npo>ni. 


|npo)p0|<(>(*"^"))i, 


( 12 ) 

UpO 

^ = 1, (13) 

ni—0 

(14) 

(15) 
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where = (1 - z) l^i)i is two-mode signal/idler emittted Hawking radiation state. The emitted 

Hawking signal/idler pairs are approximately squeezed for early time z < z*, since for long time evolution the exact 


(n) 

state in Eq.(13) does not factorize as in case of the short time state Eq.(15|. Note that notation indicates 


the probability that ni particles are emitted into the Hawking radiation signal/idler mode when there were initially 
Upo particles in the BH ‘pump’ mode. Henceforth, we shall denote \ni)i = \ni)g.j., drop the argument z on the 
probabilities, and leave implied the unoccupied vacuum signal/idler states for k' greater than the current 

timeslice considered. From Eq.([^ and Eq.([^ the state |'k(l)) is clearly normalized to unity. 

To illustrate the notation we will employ, it is instructive to write down the wavefunction at after the second 
emission event 

fT^pO fT'pO 


it,pU '■'pu '"l I 

|4'(2)) = [/p,2 |4'(1)) = ^ yP^Z"‘°\^Z"‘°^~""^\inpO-rii)-n2)p\rii)i\n2)2, 

m—O n2—0 

'/T'pO j - f^pO j - 

« \npo)p®Yl ® npo»ni,n2, ^ ^ 


(16) 


ni=0 


n2—0 




|ripo)p® 


"^pO "^1 

n2—0 


(17) 


The new feature of Eq. is that the second particle has been emitted into the second signal/idler mode with the only 
dependence upon mode 1 being that the initial number of particles in the BH ‘pump’ source is now npQ — ni, where ni is 
the number of particles that were emitted into mode 1 during the first emission event (note: ni G (1,..., N) = 1 : N). 
Again, in the short time limit Eq. 0 indicates that the emittted Hawking radiation is a succession of independent 
two-mode squeezed states in modes 1 and 2 respectively. 

Note that by utilizing a wavefunction |'k(2)) we are implicitly assuming a degree of coherency between the pump 
and the emitted Hawking radiation signal/idler modes, as exhibited in the exact states for |'I'(1)) and |4'(2)) in Eq. 
and Eq.( 16) respectively. The goal of the one-shot decoupling procedure is to decouple the emitted Hawking radiation 


modes from the pump at each emission event, while also keeping track of the finite and decreasing nature of the BH 
quantized degree of freedom np(r) that arises from the finite, though large, initial occupation number Upo 1. In 
the language of laboratory SPEC, on is making the implicit assumption that the coherency of the BH ‘pump’ source 
is shorter than the average time between emission events. We return to a discussion of this point in Section [V] 

Since each unitary emission {Up^i}i^i-,N acts for a short time At, we are continually in the short time regime z < z* 
and each emitted signal/idler Hawking radiation pair is nearly a two-mode squeezed state. However, the occupation 
number of the BH ‘pump’ mode is continually decreasing, and it is the effect of this finite nature of the ‘pump’ source 
on the total state that we wish to examine for long times (large N) as the BH evaporates. Consider the wavefunction 
|4'(7V)) after N emitted events given by the generalization of Eq.([l6| 

Kpo ripo— np0 — (ni+n2) ripo —(ni-)-...-|-niv-l) 

w«)) = E E E ■■■ 

m—O 712—0 n^—O riN—O 

N 

X |npo - («■! + • ■ • + «Af))p C) |ni)i. 


E / (n) (npo-m) (ripo-ni-ns) 

V P '^2 Pn^ ■ • ■ PriN 


N 


|«po)p o n > 


(18) 

(19) 

( 20 ) 


By construction we have ('k(iV)|'k(A)) = 1. 

Let us rewrite |'I'(fV)) as follows. We define j/ = X))n=o 7o = 0. Keeping track of the upper and lower 
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limits on each summation, we obtain the representation 

^pO "^po ^po '^po N ^ 

|vI/(7V)) = (l-z )^/2 ^ ^ ... ^ V^lnpo - J7v)^ ^ n ^ 


liz ji-l)ij 


jl=0 32=jl 03=32 3N=0 

'fT'pO 

= (1-2:)^^^ ^ V^lnpo - jAr)iV 

3JV=0 
fT-pO 

3JV=0 


i=l 

3n 3n 3n 3n N 


E E E ■•• E nI-?* 

il=0 J2=jl i3=t2 3n-i=3n-2 z = l V V ) 


(21) 

( 22 ) 


where we have defined the unnormalized state by the expression in the large square brackets in Eq.(21), and 


we have pulled the sum over the index jjv to the far left, which alters the limits of the remaining inner nested sums, 
describes the emitted Hawking radiation state with exactly jat particles (at the A^th time slice) emitted into N 


possible distinct signal/idler modes, which is in general a superposition state over many Fock states whose occupation 
numbers sum to exactly jat. 

We now wish to approximate 


for large, but finite ripo ^ 1- In Eq.(21) the factors (1 — h+i) 1/2 g^j-g 


negligibly small for all but ji ~ Upo- Even for ji = ripo the factor only contributes a (1 — « 1 + z/2 for z 1, 

the case we will consider in this work (i.e. weak two-mode squeezed states). Note the case ji = ripo corresponds to 
all UpQ particles of the BH emitted in a single burst at time Ti = i At into mode {si, ii), which constitutes a very low 
probability event for early times, but not necessarily so for longer times. Thus, to lowest order, we approximate all 
the summands by unity for all modes i = 1 : N. i.e. (1 — 2 ”!’““h-i-i)-i /2 Therefore, we have 


3n 3n 3n 


3n N 

ic’>»10 = 1: i: S'-, i: nii.-A-.).. 

jl—0 j2—jl 33—32 3 n-i—3n-2 ^ = 1 


N 


Ylbi- 

jl^j2---<jN-2<jN-l ^=1 


(23) 


with 


Jn + N-1 

3N \ 3N / \ 


(24) 


where the binomial factor in Eq.(24) counts the number states containing exactly jn Hawking radiation particles 


into N signal/idler modes, i.e. the selection oi jiq N — 1 objects taken jN at a time with repetitions. We can also 
intuitively understand the nested sum in Eq.(23) over the dummy indices ji < J 2 • ■ ■ < jN -2 < Jjv-i as the number of 


lattice points in the ‘upper diagonal’ quadrant (including the diagonal) of a iV — 1 dimension hypercube with Jn + 1 
lattices points (0, 1 ,... , jat) per dimension. 

We normalize the state as 


||.W) 

I 7Ar / 


1$ 


,{N) 


) = 


Jn + N -1 
3N 


- 1/2 


\ IM / 


. ) = 1- 

^ ripo-jN' ripQ-jNi 




(25) 


As an example 


‘^Tn=2) = (| 2 , 0 , 0 , 0 )+ 10 , 2 , 0 , 0 )+ 10 , 0 , 2 , 0 )+ 10 , 0 , 0 , 2 )+ 11 , 1 , 0 , 0 ) 

+ |i,o, 1,0) + |i,o,o, 1) + |o, 1,1,0) + |o, 1,0,1) + |o,o, 1,1)) /yio. 


JN + N -I 

Jn 


= 10, (26) 


Upon replacing jat fc (to simplify notation) as the total number of Hawking radiation particles emitted into N 
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signal/idler modes, we obtain 

^pO , - 


k I^^p0 - A:)p |$'E), 


3(^) 


fc=0 


■JT-pO - - 


= (i-zrz' 


Af k ( k + N — 1 
k 


piN) _ 
-* fc ~ 


P, 


(AT) 


Y^ripo pi,/ 

l^k'=0 ^k' 


r, (27) 


fc =0 

Af 


TipO 1 


i«po)p 0 n 


= (1 - z) E \n)i 


(28) 


(29) 


i=l 


n=0 


whe re i n the last line, the emitted Hawking radiation signal/idler field modes are in a product of N squeezed states. 
Eq.(27l with probability is one of the primary analytic result of this paper, and leads (in the next section) to 

entropy curves proposed by Page, and discussed in Blencowe and Nation [T] and Alsing [5], and examined numerically 
for this current one-shot decoupling model in Bradler and Adami Note that in the limit n^o —> oo we have 
J2T=o ~ 1 using the identity =(1~ z)~^. It is informative to compare the above state of 

the emitted Hawking radiation with the separable product of N two-mode squeezed states. 

The product of N squeezed states with = (1 — -z) (where 

\ni)i = \ni)s. \ni)ij can be written as (A^))^^ = J2T=o \^k^^) where pj.^^ "po>i Pk^^■ The density 

E OO 

k —0 2 -^ k '—0 


matrix p 


{sqzd) 




k l\'^k' I pure state contains off-diagonal y .i ^ 


as 


well as diagonal matrix elements Plf^\ but of course has only a single non-zero eigenvalue of unity, since it is a 
pure state. The bipartite state of the BH ‘pump’ mode and all signal/idler modes |d>(7V)) = X)fc=o \/^k^'^ l^po ~ 
A:)p Eq.(27l has an important difference, even in the limit n^o —>■ oo. In forming the density matrix 


Pp{N) = Tr,y[|4'(7V))(^'(Af)|] and Ps;i{N) = Trp[|^'(Af))(^'(Af)|] we pick up a factor of = ^k',k and 

p{npo — fc'jripo — k)p = 5k',k respectively, which yields the diagonal, reduced density matrices 


TlpO 


TlpO 


Ppi^) = E I’^po “ ^)p(’^po - fc|> Ps,i{N) = E 






/(AT) I 


(30) 


A;=0 


fc=0 


with non-unit probabilities P^^^ given by Eq.(27|. Thus, the underlying origin of the probabilities P^^^ arises from 


the one-shot decoupling of the sequence of separably emitted Hawking signal/idler modes, which are however, each 
separably coupled to the BH ‘pump’ mode at each time step N. This is in agreement with physical approaches to BH 
particle production advocated in previous work [3l[7l[8]. 

It is also interesting to note that the form of the probabilities Pj.^'^ given by Eq.(|27|) are reminiscent of the initially 
seeded signal states considered by both Alsing |2] and Adami and Ver Steeg |20f. That is, if instead of using the 
initial signal/idler vacuum state |0)sy = |0)s|0)j, one considers the state with Uso initial particles in the signal mode 
|'0(O)) = |npo)p |?^so)s|0)i one would obtain in the short time limit z < z* 




= (l-z) 


"sO-l-l 


k + Uso 

k 


0<z<z* 


(31) 


which is of the same form as P^’^^ in Eq.(|27|) if we take Uso = N — 1. The formal equivalence of the two probability 


distributions is the essential reason that the trilinear model considered by Alsing in [2] with only one emitted Hawking 
signal/idler mode was capable to reproducing the Page information curves for long time evolution (see Fig.14 in [2]), 
noting that the model should only be considered in the region where (inp(r)/dr < 0 holds (and whose validity stops 
once dnp(r)/dr = 0 again with r > 0). 


IV. NUMERICAL RESULTS AND FINER APPROXIMATIONS TO THE PROBABILITIES 

A. Probabilities and entropy curves 

In Figure Fig.(^ we plot plot the probabilities pI^\z) using z = 0.1 for various values of jjv = k, the collective 
number of particles emitted in all the signal/idler modes at iteration N. 
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FIG. 1: Plot of for Upo = 25 and z = 0.1. 


In Fig. 1^ we plot the von Neumann entropy = slfl\z) = — ^oSnpo+i for the reduced 

density matrix (z) = Tr^-[pp^^-] or p[^\z) = Trp[pp^^-] forripo = (5,20,50,100) (as in [3]) and TV =(!,.. ■,Nmax) 
with Nmax = (2000,4000,8000,10,000), after which we scale all graphs to N^ax = 2000. The utilization of the 
logarithms base UpQ + 1 is to ensure that the von Neumann entropies remain less than unity, for comparison. The 


S|p(NAr)] 



FIG. 2: Sp^\z) vs N for Upo = (5,25,50,100) and z = 0.1. Entropies are computed with 

entropy curves in Fig.([^ have the characteristic features of (i) an initial value of zero, appropriate for an initial 
product state at = 0, (ii) a rise to a peak value less than unity, and (iii) a tapering off towards zero for long-times 
(large N). 

An examination of the probabilities for various values of N for Upo = 25 Fig.(|^ reveals that they are peaked 
near fc = 0 for early times, when few Hawking particles are in the signal/idler modes and the state is essentially a 
separable product of squeezed states with the BH ‘pump’ mode occupation number near its initial value of ripo- As 
time, i.e. N, evolves this peak moves steadily from low values of k to high values of k, and for long times clusters 
about k = Upo, where the BH has essentially evaporated. In this longtime regime, we would expect BH ‘pump’ to 

be in the state jripo — fc « 0)p with the emitted signal/idler held approximately in the state ^), i.e. again, an 

approximate product state for |4'(A^)). As such, we would expect the entropy curves for longtimes to again be zero. 
We address this in the next section by further approximating the probabilities ■ 


B. Finer approximations to the probabilities 


We desire to further approximate the the unnormalized state 14*^^) dehned by the expression in the large square 


brackets in Eq.(21). In particular, in the previous section we had approximated the terms (1 — 2 ;"po-ii-i+i)-i /2 i^y 
unity. Recalling that jo = 0, we can factor out from all the nested sum an overall constant term (1 — 2;"po+i)“i/2 —1 
for z 1 and any reasonable sized value of UpQ. The remaining factors have to be summed from j^ = ji_i : = fc. 
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FIG. 3: vs k for ripo = 25, 2 = 0.1 and N = (10,100,200,300,500,1000). 


in succession from the inner summations, outwards. These complicated nested sums are what led to the numerical 
lattice-path approach of Bradler and Adami Here, we make the simplified, but reasonable approximation that 
(1 — is dominated by its largest contribution ji = k from the upper limit of the summation, yielding 

(1 — which can then be factored out of all the nested summations, except the outermost one over k 

itstelf. Since there are iV — 1 summations at time N we obtain the slightly refined approximation to the probabilities 


pW 


( 1 - 2 ) 


N 


(1 — 


k + N -I 
k 


.'W 


p(N) 

_ 

p{N) ■ 

2^k'=0 ^k' 


(32) 


Eq. (321 constitutes the second primary analytical result of this present work. 


The effect of the extra factor (1 — 


S[p(NAr)] S[p(NAt)] 
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— 25 

— 50 

— 100 


FIG. 4: Entropies with the extra term (1 - 2 "po-fc+i)min{jv-i, 
(right) ripo = 25. Entropies are computed with 


in Eq.(32l with various values of imax. for (left) ripo = 10, 


2 "po k+i'^N 1 piN) replacing it with (1 — i.Jmax i) varying the value of 

imax < N- The value of imax sets how many terms (1 — z"po-fi-i+i)-i/2 TV — 1 nested sums in Eq.(|2l[) that 

we do not approximate as unity. This is shown in Fig.Q for the cases of Upo = 10 (left) and Upo = 25 (right jTThese 

figures show how the additional factors of (1 — in Pjf^^ brings down the tail of entropy distribution S to 

zero for longtimes, while leaving the short time (small N) portion of S essentially unaltered. 
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C. Page Information 


To investigate how the information emerges from the BH as it evaporates, the Page information^ we follow Page’s 
1993 paper [5] (and [HE]) and define the information I as 

I{t) = SthermaliN) - (33) 

Here Sthermai is the effective thermal distribution pthermaiizthermai) with probability distribution given by = 

(1 - Zt/iermai) Eq. ([^ , with Zthermai = ^^^/(hsy + 1), and = X)fc=o ^ Computed from Psg,{T) = 

Tvp [|\1'(A^))(4'(A^)|]. Lastly, from Eq.([^, the (squeezing) rapidity z is dehned as z = tanh^(T) with t = ^Upo r t. Tak¬ 
ing tff = NAt = N/Njnax for N = 1 : Nmax we utilize (with r = 1) z = tanh^ [tanh“^(zma£!;) tanh(y/rqlo 
with 0 < Zmax < 1 where < Nmax determines a controllable build-up time from z = 0 to z = Zmax- Since 

each unitary {Up^i\\i-i-,N acts over a small time, we allow the rapidity to gradually build up from z = 0 : Zmax 1, 
assuming weak squeezing in the emitted Hawking signal/idler radiation pairs by the BH. (Note: similar curves are 
obtained by simply setting z = Zmax for all N). 



FIG. 5: Plots of entropy S (black, solid), effective Sthermai (gray, solid). Page Information I (black, dashed) and total number of 
emitted Hawking particles in all signal/idler modes n^ t/npo vs time N for ripo = 25 and Zmax = 0.1, = 200, imax = 50), 

using probabilities (left) Pj^^^ in Eq.|27L and (right) in Eq.(32l. Entropies are computed with log^^^^^j^. 


In Fig.([^ we show plots of the entropy S (black, solid), effective Sthermai (gray, solid), Page Information I = 
Sthermai — S (black, dashed) and total number of emitted Hawking particles in signal/idler modes fig ifnpo vs time 

N for Upo = 25 using probabilities (left) Pj^^^ in Eq.(27l, and (right) Pj^'^ in Eq.(32). Entropies are computed with 
logripo+i that all graphs have maximum value of unity, for comparison. Both curves show that for early times 
(small N) S « Sthermai SO that the Page information / is flat with, with very small slope. As time progresses, / 
begins to grow, as rapidly increases, and the BH begins to evaporate. For h^y/ripo >1/2 there are less particles 
in the BH ‘pump’ mode than have been emitted into all the Hawking radiation signal/idler modes and S begins to 
decrease. Using pjf^'^ from Eq.(27) S in Fig. ([5|-( left) decreases very slowly as |4'(Af)) approache^a nearly separable 
state Iripo ~ k ; 


0)pl<J>fc 


(TV) 


Tipo'^s.iy and the information slowly saturates to a value of unity. In Fig.(|5|)-(right) the effect 


of the refined probabilities Pj'^'^ from Eq.(32) reduce S more rapidly to zero, appropriate for a final separable state, 
with a commensurate faster saturation of the information to unity. 

In Fig.([^ we show the similar entropy and information curves for Upo = 100, where the initial flatness of the 
information / is more pronounced. Fig.([^-(right) and Fig.Q constitutes the main results of this present work. Note 

f k + N — 1 

that a brute force summation of all the terms in Eq.(|21|) would involve the addition of on the order of ( ^ 


summands, which equates to 10"^^ and 10^®^ terms for k = Upo = 25, N = 500 and k = ripo = 100, N = 2500 for 
Fig.® -(right) and Fig.([^ respectively, which is impractical. While most of the summands would be negligibly small 
to warrant approximating to zero, a reasonable estimate of only fc = 10 nonzero terms per sum would still lead to the 
prohibitive total number of nonzero summands of 10^° and 10^^ for N = 500 and N = 2500, respectively. Hence, the 
necessity for the analytic approximations to the probabilities given by Eq.(27) and Eq.(32). 

As explored in Alsing [2], the intial BH ‘pump’ might also be represented by a coherent state |a) = 
g-|“l /2 a"|n)/v/n! such that ap\a) = q;|q;), the quantum state that most approximates a classical state with 
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n=100 

^thermal » n/n[p(NAT)] 



^thermal 

s 


“““ 1 ^thermal ^ 

-n/n 


FIG. 6: Plots of entropy S (black, solid), effective Sthermai (gray, solid). Page Information I (black, dashed) and total number 
of emitted Hawking particles in signal/idler modes n^ -ifnpo vs time N for Upo = 100 and Zmax = 0.1, = 10^, imax = 50), 

using probabilities in Eq.(32l. Entropies are computed with log„^^_,_j^. 


mean particle number = UpQ. This leads to the state |tE'(A^)) « X]^o=o \/Sfc=o l^po — 

^)p probabilities [52] which tends to smooth out the entropy curve S as 

shown in Fig. ([^-( right) by averaging original probabilities over the initial coherent state probability distribu¬ 
tion ■ The main point is that is well approximated by with = UpQ. 



FIG. 7: (left) Initial coherent state probabilities for = 25, (right) entropy curves S for probabilities (black, solid) 


in 


Eg. ([ 2 ^, (gray, solid) in Eg.(32 1 , and (gray, dashed) (q) Entropies are computed with In. 


V. DISCUSSION 

Here we propose that the process of BH evaporation has a very strong analogy to the process of spontaneous para¬ 
metric down conversion (SPDC) [211101 with the emitted Hawking radiation acting as the spontaneously signal/idler 
pairs generated by the trilinear Hamiltonian Eq.([^ driven by the BH modeled as a depleted pump laser source. In 
Fig.([^(a) we illustrate a laser excitation source powered by an unlimited power source (e.g. A/C wall socket) driving 
an optical storage cavity. This optical cavity, with non-zero mirror transmission subsequently pumps a non-linear 
crystal in which a pump photon of frequency uip spontaneously down-converts two photons from the vacuum, the 
signal at frequency uig and the idler at frequency such that energy is conserved, ujp = uis Puji- As long as the power 
source is unlimited, signal/idler pairs are continually created as squeezed state signal/idler pairs, and the pump can be 
treated as quantized, though with large enough occupation number UpQ to essentially be treated as a (non-depleting) 
constant. This is the analogy of a non-evaporating BH generating purely thermal Hawking radiation. The laser 
excitation source filling the storage cavity can loosely be thought of as the analogy of the BH formation, while the 
storage cavity can be thought of as some region of finite width about the BH horizon (see Discussion in [2]). 

Fig.(|^(b) shows essentially the same set up as in Fig.([^(a) except for one crucial difference: the initial laser 
excitation source is driven by a power source of limited energy, here illustrated as a battery. Again the BH formation 
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FIG. 8: Analogy of BH evaporation to SPDC process. 


fills up the storage cavity, but now with a large but finite number ripo of quantized particles which can still drive the 
SPDC Hawking radiation production process. In this later case, the storage cavity will eventually deplete itself (e.g. 
finite battery life). The signal/idler pairs will still be produced as essentially squeezed states at each SPDC excitation, 
with each emitted pair entangled with the BH ‘pump’ source at the given time step N (but not at subsequent time 
steps). The analogy of the non-linear crystal responsible for the SPDC generation of Hawking signal/idler pairs is the 
curvature distortion near the BH horizon, as discussed in EHH]. 



N = 6 N = Z /V = i_ 

^- *£•, ^ 


(a) Long-Coherence-Length Pump Photons - Single Frequency Laser 


Cn 



^ ' ^N=l 

^N=2 



{b) Medium-Coherence-Length Pump Photons - Quasi-Continuous Wave Laser 


Cn 





CO 


Short-Coherence-Length Pump Photons - Short Pulse laser 


FIG. 9: Gonsideration of coherence length of BH ‘pump’ source particles. 


By simply writing down a pure state wavefunction |f['(7V)) in Eq.(22| to describe our quantum state, one has 


implicitly assumed some level of coherency of the pump source which in turns effects the level of entanglement of 
the emitted signal/idler particles with the BH ‘pump’ particles, which without loss of generality, and for comparison 
with our laser SPDC analogy, we treat in this discussion as photons. This is illustrated in Fig.(|^. In Fig.Q(a) 
we illustrate a pump laser source of infinite coherence length ~ cjALOp. Here, Awp is the bandwidth of BH 
‘pump’ photons which for Fig.® (a) is considered very narrow, Acop ^ ojp. Consequently, all signal/idler SPDC 
pairs generated are coherent with the pump simultaneously, and therefore, indirectly with each other, causing various 
degrees of entanglement as explored by Alsing [2]. This is the case considered by Nation and Blencowe [1] and Alsing 
[5], which can also be equivalently considered as having the non-linear material/curvature distortion occurring within 
the storage cavity. As the coherence length of the BH pump photons decreases (Awp increases) Fig.([^(b), less and less 
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SPDC generated pairs are coherent with with the pump, and hence indirectly with each other. If the coherence length 
of the pump photons is less than the average time Ar {AN = 1 in our model) for SPDC signal/idler pair generation, 
as illustrated in Fig.([^(c), each generated SPDC signal/idler pair is only coherent with the BH ‘pump’ until the 
next pair is generated. This is the one-shot decoupling model of Bradler and Adami [3], and explored analytically 
in this present work. As discussed earlier, this model is commensurate with physical models for Hawking radiation 
production iSHg in which the curvature distortion at the BH horizon (of shrinking radius) creates the signal/idler 
pairs, which then propagate away from the region of generation, since the BH horizon decreases, shifting the region 
of curvature distortion for subsequent pair production. 

As evidenced by the finer approximation utilized in the previous section to produce the probabilities in Eq. (321, 

the terms (1 — in the temporally nested sum in Eq.(22) describe the detailed probability 

structure of the total BH/Hawking radiation state |'I'(7V)). Hence, further refinements than simply considering their 
maximum contribution from ji = k are warranted, yet difficult to implement analytically due to the large number of 
temporally ordered nested sums. Thus, it appears that some blend of an analytical approach detailed here, coupled 
with a numerical approach, as advocated by the lattice path method Bradler and Adami [3] might prove fruitful in 
gaining further physical insight the nature of the one-shot decoupling bipartite state |'I'(A)), and subsequently the 
detailed nature of the Page information of the Hawking radiation from an evaporating BH. 

Considering extensions of this present work, we note that our model utilized a monochromatic, single frequency BH 
‘pump’ source to drive the signal/idler generation i.e. a delta function frequency distribution for the pump. A more 
general model, which might encapsulate a more realistic physical scenario, would be to model the BH ‘pump’ source 
not as a single frequency, but as a collection of frequencies over some bandwidth, ie a pulse source. Our model also 
described the emitted Hawking radiation signal/idler pairs as single frequency modes, which is analogous to a long 
non-linear crystal with strict phase matching conditions mm- A more general model could incorporate a frequency 
bandwidth for the emitted Hawking radiation signal/idler pairs, similar to a shorter non-linear crystal with a less 
restrictive phase matching condition. 
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